Electronic Orbital Currents and Polarization in Mott Insulators 
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The standard view is that at low energies Mott insulators exhibit only magnetic properties while charge de- 
grees of freedom are frozen out as the electrons become localized by a strong Coulomb repulsion. We demon- 
strate that this is in general not true: for certain spin textures spontaneous circular electric currents or nonuni- 
form charge distribution exist in the ground state of Mott insulators. In addition, low-energy "magnetic" states 
contribute comparably to the dielectric and magnetic functions eik{uj) and ^ik[uo) leading to interesting phe- 
nomena such as rotation the electric field polarization and resonances which may be common for both functions 
producing a negative refraction index in a window of frequencies. 



Mott insulators are the paradigm of strongly correlated ma- 
terials. Their minimal Hamiltonian is the Hubbard model 
which includes a hopping term, t, and on-site Coulomb in- 
teraction U . At half filling (one electron per site) and in the 
large U jt limit, each site is occupied by a single electron 
to avoid the strong on-site repulsion. The charge becomes 
localized by this mechanism and the low energy properties 
are described by the remaining spin degrees of freedom. For 
this reason, Mott insulators with large U jt have been tradi- 
tionally considered as materials which have only magnetic 
properties at low energies due to their spin moments. De- 
spite this common conviction, we will show here that certain 
ground states of Mott insulators exhibit real electric currents 
in loops (orbital currents) that produce orbital magnetic mo- 
ments, while others show modulation of electron charge (po- 
larization). Consequently, spins in Mott insulators are cou- 
pled not only to dc magnetic fields, but also to dc electric 
fields, and it is possible to have magnetically driven electronic 
ferroelectricit y H-i^ii3.i4 Moreover, nonvanishing matrix ele- 
ments of the polarization between the ground state and excited 
magnetic states result in a nonvanishing contribution to e^fe [lo) 
at low energies with optical strengths comparable to those of 
[lik (t^) for the diagonal and off-diagonal elements. Therefore, 
rotation of the electric field polarization is a characteristic sig- 
nature of spin textures with orbital currents. 

The apparent contradiction between the insulating nature 
of the system and the existence of non-zero orbital currents 
is resolved when we notice that electrons are not completely 
localized on their ions for finite U jt. In fact, the effective 
Heisenberg interaction, J o^i^ /U, results from a partial delo- 
calization: an electron gains kinetic energy by "visiting" vir- 
tually a neighbouring site, but this only occurs if the spins are 
opposite on both sites (Pauli principle). Simlarly, the elec- 
tron moves along a closed loop generating local currents that 
depend on the spin structure along the loop. 

We start by considering a half filled Hubbard model on a 



general lattice: 



where sites are labeled by indices {ci^) is the cre- 

ation (annihilation) operator of an electron with spin a on 
site i and rii = J2a Aa^ia '^he number operator. The low 
energy spectrum of this model is described by an effective 
Heisenberg spin Hamiltonian, H, which is obtained by the 
usual degenerate perturbation theory mtjU << 1. H acts 
on the low-energy subspace where all the sites are singly oc- 
cupied. Consequently, it is expressed in terms of the spin op- 
erators: Si' = J^ij. V AiiP^liv'^'iv^ where g"^ are the Pauli ma- 
trices and r\ ~ {x, y, z}. The expression of H to order is: 
^^'^ - (SrS,- 1/4), with = itl^/U. In general, 

any physical operator, O, has an expression in terms of spin 
operators, O, that results from the application of degenerate 
parturbation theory. 

We will consider first the current operator 



hr. 
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between sites i and j. Since the shortest loop is a triangle, 
the lowest order finite contribution to the current operator is 
/U"^ and contains the product of 3 spin operators. The cur- 
rent is a scalar under spin rotations and it is odd under time 
reversal and under spatial inversion. The only possible ex- 
pression involving three spin operators is the so-called "scalar 
spin chirality" operator. Using perturbation theory'-, we find 
that the contribution to the current in the bond 1 , 2 from the 
triangle 1-2-3 is: 



ri2 24eii2t23t3i r„ ^ „ , 

-1-12,3 = T -rp\ 31 X 532 • 
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The quantity Xy.fe = [Si x Sj] • Sfc called scalar spin chi- 
rality was introduced previously in numerous discussions of 
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magnetic systems." It was invoked in the theory of anyon 
superconductivity^ and for describing properties of triangu- 
lar and kagome magnets^. Scalar chirality produces a novel 
Berry-phase mechanism for anomalous Hall effect^. Chiral- 
ity can lead to new universality classes in phase transitions^ 
and special "chiral glass" phases in disordered systems^. It 
can also modify quantum tunneling in magnetic molecules^ 
and it appears in the electronic contribution to the Raman 
scatteringiS. Despite this broad interest, the physical meaning 
of scalar chirality remained unclear. Eq.© shows that scalar 
spin chirality corresponds to orbital currents running in the 
low energy states of Mott insulators. These currents produce 
orbital magnetic moments hijk oc Xij,kZ, where z is normal 
to the plane of the triangle. We note that orbital currents only 
appear for noncoplanar spin structures, as is proportional 
to a solid angle formed by the spin vectors Si, S2 and S3. 

Normally, orbital magnetism leads to a paramagnetic re- 
sponse. However, since the coupling of magnetic field to or- 
bital moments is weak (oc /U"^) compared to the spin Zee- 
man coupling, the dependence of on the external mag- 
netic field is mainly due to changes in the spin configuration. 
For this reason can increase or decrease with field. For 
instance, the coplanar 120° spin structure in the easy-plane 
triangular lattice is turned into an "umbrella" pattern when 
a magnetic field perpendicular to the easy-plane is applied, 
see Fig. [T^. The orbital moment on a triangle ijk becomes 
nonzero and its absolute value first increases with field but fi- 
nally decreases to zero when the soins become fullv oolarized. 




FIG. 1: (a) Ground states with nonzero electric current of the C3 
invariant Heisenberg triangle. The circular arrows indicate circular 
currents, and (+, — ) the sign of a scalar chirality. (b) Examples of 
magnetic states with nonzero polarization. The two spins inside the 
oval form a singlet state. The unpaired spin can be up or down 

The orientation of the orbital moments depends on the signs 
of tij . For lattices, the current on a given bond is the sum of 
the loop currents in all the triangles to which that bond be- 
longs: ly ~ i^ij/Tij) J2k ^ij k - In some particular cases, the 
net current of the bond is zero because different contributions 
cancel each other. In this situation a net current exists on the 
surface of two-dimensional lattices. 

For instance, for regular spin tetrahedra (building blocks 
of pyrochlore or spinel lattices) with strong uniaxial 
[±1, ±1, ±1] anisotropy, the structures "four in" or "four out" 
do not have net currents, see Fig|2^. However, the structures 
"two in - two out" (spin ice), or "three in - one out", which 



could be stabilized by an external field, have nonzero net or- 
bital currents as shown in Figs|2j3,c. Similarly, two widely 
discussed structures in kagome lattices are those with a homo- 
geneous vector chirality Si x S2 + S2 x S3 + S3 x Si (q = 
structure), and with staggered vector chirality (V3 x \/3 
structure) ^''^'^^ . In both cases, there is easy plane anisotropy 
and the umbrella structure induced by a magnetic field per- 
pendicular to the lattice has a nonzero orbital moment. As 
shown in Fig 13 the pattern of currents and orbital moments 
is uniform in the first case, and staggered for the latter case 
(despite the fact that the net spin moment is the same). The 
coupling of a net orbital moment to an external magnetic field 
favors the uniform state. 

In a similar way we derive an expression for the projected 
local electron number operator fii. This operator is a scalar 
under rotations in spin-space, i.e., it must be a function of 
the combinations S^ • Sj . The first non-zero contribution to a 
deviation Shi from unity is 

Sfi, = ni - 1 = gMp!^ [Si . (S2 + S3) - 2S2 • S3] . 

(4) 




FIG. 2: Non-coplanar spin configuration in a pryrochlore lattice, 
a) Four spins point inwards along the principal diagonals; the net 
currents are zero, b) Three spins point inwards while the other one 
points outwards leading to a net current circulating in the opposite 
triangle, c) Two spins point inwards and the other two point out- 
wards. The orbital current circulates in a loop formed by four edges 
of the tetrahedron. 

A similar expression holds for the charges at sites 2 and 
3 after a cyclic permutation of indices. The spin structure 
of the charge operator is uniquely fixed by the invariance of 
hi under the time reversal symmetry and the interchange of 
sites 2 and 3, as well as by the conservation of total charge 
of the triangle: Y^^=i — The redistribution of charges 
within triangles induces electric dipoles in magnetic states: 
Pe = e ViShi. These dipoles can appear spontaneously or 
can be induced by a magnetic field, see FiglTJ). In particular, 
a triangle with the classical coplanar 120° ordering of spins 
does not have charge redistribution in the ground state. How- 
ever, for an easy-plane anisotropy, an in-plane magnetic field 
perpendicular to the bond 2-3 cants the spins in such a way 
that (Si • S2) = (Si • S3) becomes larger than (S2 • S3) induc- 
ing an electric dipole moment in the field direction. Thus, a 
nonzero electric polarization of pure electronic origin appears 
in such nonuniform spin configuration. 

The projected dipole and current operators are identically 
zero for bipartite lattices with nearest-neighbor hoppings— . 
This results from the invariance of H under the product of 
particle-hole and t ^ —t transformations, while P and I are 
odd under this transformation. 
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electric and magnetic fields. The matrix elements of P be- 
tween the ground state |0) and excited magnetic states jri) de- 
fine the contribution of these states to the dielectric function. 



FIG. 3: Spin configurations on a kagome lattice that lead to a current 
or charge ordering in an applied magnetic field, a) The "umbrella" 
phase induced by a field perpendicular to the plane has a uniform cur- 
rent ordering. The ± signs denote directions of current, b) The same 
as a) for staggered current ordering, c) Spin ordering induced by 
field for a Heisenberg model on a kagome lattic o'^i'^ . The elongated 
ovals indicate a resonant valence bond state on the corresponding 
hexagons, while the little circles represent spins that are polarized 
along the field direction. This structure is accompained by charge 
ordering: the charges on sites belonging to hexagons and those on 
isolated sites should be different. 



We will discuss now some consequences of the obtained 
results. The electric dipole induced by virtual electron hop- 
ping has the same form as the one resulting from the de- 
pendence of the exchange constants on ion displacements 
in a magnetically ordered state (magnetostriction): w 
Jij{0) + u„ • V„Jy . Minimizing the sum of the magnetic 
energy ^ JijSi ■ Sj and the lattice distortion energy respect 
to Ui, we find that the resulting electric dipole of a triangle 
is expressed in terms of scalar products • S ^ . Due to the 
symmetry considerations discussed above, the spin structure 
of the dipole is the same as in Eq.(|4l), while the coefficient is 
~ e |V J| /K, where K is the lattice spring constant. Elec- 
tronic dipoles, Pe, corresponding to Eq. (HJi, together with 
the spin-dependent dipoles originating from magnetostriction 
lead to the coupling between spins and electric field (eZ; is 
the charge of an ion /): 

= -P • E, P = Pe + ^ eZiKf\S, ■ - 1/4). (5) 

ijl 

A simple example of a spin-driven charge modulation 
is given by the the 1/3 plateau phase of the S = 1/2 
kagome lattice. This state has the local structure shown in 
Fig IS, with a resonating singlet state on the hexagons and 
up-spins in between^^, and could be a long-range ordered 
valence-bond crystal^Z.. A similar situation arises forthe sev- 
eral magnetization plateaus in the Shastry-Sutherland system 
SrCu2(B03)2^. The states at each plateau consist of or- 
dered arrays of singlet and triplet dimers which according to 
(IHi should lead to a spin driven charge density wave. There 
are also systems consisting of isolated triangles with long- 
range magnetic ordering, e.g. La4Cu3MoOi2^°. According to 
Eq. (lUi, the magnetic structure found in Ref. 20 should have 
nonzero total electric polarization (multiferroic behavior). 

Next we consider the response of the Mott insulator to ac 



eik{uj) = eoOik + -y- 



E 



uJno{0\P,\n){n\Pk\0) 



(6) 



at r = and at frequencies well below the frequencies of op- 
tical phonons. Here S ^ 0, huj„Q = En — Eo,H\n)^ En\n). 
Further, eo is the contribution of all the other high frequency 
modes, V is the total volume. The expression for the mag- 
netic response function, ^ikiuj), is obtained by replacing P 
with guB^ (we neglect the effect of L relative to the spin 
contribution and the difference between S and S). 

In absence of spin-orbit coupling, the states \n) have well 
defined total spin S and z-projection 5*2 (they are eigenstates 
of the isotropic Heisenberg Hamiltonian). The operator P pre- 
serves these quantum numbers. In contrast, and 5*1, connect 
states with different total spin. Therefore, the excited states 
that contribute to eik{uj) and Hik{'jj) are in general different 
and their resonances are different too. However, in the pres- 
ence of spin-orbit coupling some resonances may be common 
for both as we will see below. It is important to note that the 
matrix elements of P are of order 8eat^/U^, i.e., about the 
same order of magnitude as the matrix elements of .g/isS for 
J ~ 100 K. Here a is a characteristic interatomic distance. 
Hence, the response of a Mott insulator to an ac electric field 
may be similar in magnitude to the response to an ac magnetic 
field. We note that if |0) is an eigenstate of with nonzero 
eigenvalue (orbital currents), the matrix elements (0|Pa;|0) 
and (0|Pj,|0) are simultaneously nonzero leading to circular 
dichroism or rotation of the electric field polarization. This 
rotation is almost the same as Faraday rotation induced by 
spins on the ac magnetic field polarization. Hence, detecting 
of orbital currents ("scalar spin chirality") is possible by mea- 
suring the rotation of the electric field polarizationi.22, To see 
that {0\Pi\n) 7^ for i = x,y if there is a net orbital current 
in the state |0), we note that for electrons moving on a ring 
(or loop in general) the orbital moment is the conjugate 
variable of the angle ip and Px oc cos Lp, while Py oc sin ip. 

The single equilateral triangle of = 1/2 spins with an- 
tiferromagnetic Heisenberg and spin-orbit interactions pro- 
vides the simplest realization of our results. There are many 
compounds, known as trinuclear spin complexes, that contain 
such isolated triangles. Triangular clusters exist in magnetic 
molecules like "V15" K6[Vl^As6042(H20)] • 8H20i^ or 
form well-ordered solidsaS. In the absence of Dzyaloshinsky- 
Moriya (DM) coupling, the ground state is a quartet with total 
spin = 1/2. The higher energy S = 3/2 quartet is sepa- 
rated from the ground state quartet by the gap 3 J/2. The four 
lowest degenerate states, |x, cr) , are labeled by "spin chirality" 
X = ±1 (which plays the role of the z-projection of a pseu- 
dospin variable T) and spin projection of Sz {cr =t or D- The 
full space of the quartet can be presented as a direct product of 
spin and pseudospin subspaces. It turns out that P^. cx ST^ = 
2S2 • S3 - Si • (S2 + S3), Py OC VSTy = Si(S2 - S3) and 
and L2 oc ^/3Tz/2 — Xi2.3^ where the Pauli matrices T oper- 
ate in the pseudospin subspace and obey SU(2) commutation 
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relations. This is a consequence of the fact that current and 
polarization are associated with conjugate variables (angular 
momentum and phase). Thus matrix elements of both and 
Py are nonzero between states with opposite orbital momenta 
and spin chiralities. Note that the eigenstates of break time 
reversal symmetry (see Figll^), while eigenstates of Px ± iPy 
break the spatial C3 symmetry (see FigHJ)). 




FIG. 4: Energy levels of spins on triangles with the account of ex- 
change, Dzyaloshinsky-Moriya and Zeeman interactions as a func- 
tion of applied dc magnetic field. Blue (green) arrows show transi- 
tions induced by ac electric (magnetic, EPR) field. 

In some real systems like V15 (see Ref.[l^ the lattice sym- 
metry allows for a nonzero DM coupling i?DM = '^ij ' 
[Si X Sj] . The terms that mix the 5 = 1 /2 and S* = 3 /2 states 
(with in-plane components of the vector D^) are relatively 
small. On the other hand, the term plays the role of a spin- 
orbit coupling between the spin and the orbital moment 
and splits the ground state quartet into two doublets {|+, |), 
|-, J,)} and {| + , I), |-, t)}, separated by an energy A = Dj. 
Consequently, the system exhibits the following properties: I. 
In absence of a static magnetic field, the electron paramag- 
netic resonance (EPR) spectrum exhibits the same resonance 
frequency, ujq = A/fi,, as the dipole-allowed microwave ab- 
sorption with similar intensities. In a static magnetic field, 
B, the absorption frequency due to the ac electric field re- 
mains the same while the EPR frequency splits linearly in B, 



see FigU II. Slightly below uiq, the contribution of this res- 
onance to both eii{uj) and i-iuiu;) for i ^ x,y is negative and 
both will be negative if dissipation is weak (negative refrac- 
tion index)22. III. Off-diagonal elements exy{^j->) and eyx{oj) 
are nonzero at low temperatures in the presence of the mag- 
netic field which splits lowest doublet resulting in a strong 
rotation of the electric field polarization at frequencies of or- 
der ujQ. IV. The electric field causes transitions between states 
with the same total spin. Therefore, eih{uj) changes with the 
ground state magnetization until the contribution of magnetic 
states vanishes when all spins become aligned in the field di- 
rection. Hence, measurements of the dielectric function, in- 
cluding the static case, provide information about the structure 
of magnetic spectrum. 

While the states with nonzero orbital currents or polariza- 
tion are degenerate in a single Cs-invariant triangle, this may 
not be the general case of infinitely large systems containing 
triangles in their structures. Different symmetries can be bro- 
ken in these systems. If the resulting spin ordered state is such 
that (xij.k) 7^ 0, Eq. (O implies that the spin ordering is ac- 
companied by an ordering of orbital moments. 

In conclusion, magnetic states of Mott insulators show elec- 
trical properties which distinguish them from standard band 
insulators. These states can exhibit electric dipole moments, 
which gives a purely electronic mechanism of multiferroic be- 
havior. Spin states also contribute to the low-frequency op- 
tical properties, such as absorption by magnetic excitations 
well below the gap for single-electron excitations (Hubbard 
gap). An even more striking property of Mott insulators is the 
presence of orbital electronic currents and the corresponding 
orbital moments which can be detected by measuring the re- 
sulting rotation of the electric field polarization or by nuclear 
magnetic resonance. 

The authors thank A. Saxena for useful discussion. LANL 
is supported by US DOE under Contract No. W-7405-ENG- 
36. The work of D.Kh. was supported by the DFG via SFB 
608 and the European project COMEPHS. 



' Also at the Loughborough University, Loughborough, UK. 

^ X.G.Wen, F.Wilczek and A Zee, Phys. Rev. B 39, 1 1413 (1989). 

^ Actually Xi2,3 characterizes the magnitude of an orbital current 

as shown in the paper and this object is achiral according to the 

strict mathematical definition of chirality. 
" T. Momoi, K.Kubo and K.Niki, Phys. Rev. Lett. 79, 2081 (1997). 
^ W. Schweika, M.Valldor and PLemmens, Phys. Rev. Lett. 98, 

067201 (2007). 

* Y.Taguchi et al.. Science 291, 2573 (2001) and references therein. 
' H. Kawamura, J. Phys. Condens. Matter 10, 4707 (1998). 

* H.Kawamura, Phys. Rev. Lett. 68, 3785 (1992) 

' H.Nojiri, E.Ishikawa and T.Yamase, Progr. Theor Phys. Suppl. 
No. 159, 252 (2005) 

B.S.Shastry and B.I.Shraiman, Phys. Rev. Lett. 65, 1068 (1990); 

PE.Sulevski et al., Phys. Rev. Lett. 67, 3864 (1991). 
" Y.Tokura, Science 312, 1481 (2006) 

S.-W.Cheong and M.V.Mostovoy, Nature Materials 6, 13 (2007). 
" D.I.Khomskii, Journ. Magn. Magn. Mater. 306, 1 (2006). 



C. D. Batista, Phys. Rev. Lett. 89, 166403 (2002). 
" Bogoliubov came very close to this result long ago, by using a 
perturbative approach with respect to the small parameter t/U to 
describe the correlated metallic state. N. N. Bogolyubov, Lectures 
on Quantum Statistics, Gordon and Breach, New York, p. 120 
(1967). 

M. E. Zhitomirsky and H. Tsunetsugu, Progr. Theor Phys. Suppl. 

No. 160, p.361 (2005). 
" D. Cabra et al, Phys. Rev. B 71, 144420 (2005). 

L.N. Bulaevskii, JETP 24, 154 (1967). 
" A. Tarantul, B. Tsukerblat, and A. Muller, J. Chem. Phys. 125, 

054714 (2006). 

YQiu et al., Phys. Rev. B 71, 214439 (2005) 
^' S. Miyahara and K. Ueda, J. Phys.: Condens. Matter 15, R327- 
R366 (2003). 

CM. Varma, Phys. Rev. Lett., 61, R3804 (2000). 
Veselago V.G. Soviet Physics Uspekhi, 10, 509 (1968). 



